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ABSTRACT 

We have performed 2.5-dimensional general relativistic magnetohydrody- 
namic (MHD) simulations of the gravitational collapse of a magnetized rotating 
massive star as a model of gamma ray bursts (GRBs). The current calculation 
focuses on general relativistic MHD with simplified microphysics (we ignore neu- 
trino cooling, physical equation of state, and photodisintegration). Initially, we 
assume that the core collapse has failed in this star. A few M black hole is 
inserted by hand into the calculation. The simulations presented in the paper 
follow the accretion of gas into a black hole that is assumed to have formed before 
the calculation begins. The simulation results show the formation of a disk-like 
structure and the generation of a jetlike outflow inside the shock wave launched 
at the core bounce. We have found that the jet is accelerated by the magnetic 
pressure and the centrifugal force and is collimated by the pinching force of the 
toroidal magnetic field amplified by the rotation and the effect of geometry of 
the poloidal magnetic field. The maximum velocity of the jet is mildly rela- 
tivistic (~0.3c). The velocity of the jet becomes larger as the initial rotational 
velocity of stellar matter gets faster. On the other hand, the dependence on the 
initial magnetic field strength is a bit more complicated: the velocity of the jet 
increases with the initial field strength in the weak field regime, then is satu- 
rated at some intermediate field strength, and decreases beyond the critical field 
strength. These results are related to the stored magnetic energy determined by 
the balance between the propagation time of the Alfven wave and the rotation 
time of the disk (or twisting time). 
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1. Introduction 

Gamma ray bursts (GRBs) are one of the most enigmatic and most energetic events in 
the universe (e.g., Klebesadel, Strong, & Olson 1973; Fishman & Meegan 1995; van Paradijs, 
Kouveliotou, & Wijers 2000). GRBs and the afterglows are well described by the fireball 
model (e.g., Piran 1999; Meszaros 2002), in which a relativistic outflow is generated from a 
compact central engine. Rapid temporal decay of several afterglows is consistent with the 
evolution of a highly relativistic jet with bulk Lorentz factors ~ 10 2 — 10 3 (e.g., Sari, Piran, 
& Halphen 1999). The formation of relativistic jets from a compact central engine remains 
one of the major unsolved problems in GRB models. 

What is the central engine of GRBs? From recent observations, some evidence was 
found for a connection between GRBs and the death of massive stars. Analyses of host 
galaxies show a correlation between star-forming regions and the position of GRBs inside 
the host galaxy (Bloom, Kulkarni, & Djorgovski 2002). A "bump" resembling that seen in 
light curves of Type Ic supernovae has also been detected in the optical afterglow of several 
GRBs (e.g., Bloom et al. 2002; Garnavich et al. 2003). GRB 980425 has been associated 
with an optical supernova, SN 1998bw (e.g., Iwamoto et al. 1998; Woosley, Eastman, & 
Schmidt 1999). From the early-phase observation of a very bright afterglow (e.g., Uemura et 
al 2003; Price et al. 2003) it has been revealed that GRB 030329 is likely to be in association 
with SN 2003dh (Stanek et al. 2003; Hjorth et al. 2003). Several authors (Frail et al. 2001; 
Panaitescu & Kumar 2001; Bloom, Frail, & Kulkarni 2003) have studied beaming angles 
and energies of a number of GRBs. They have found that central engines of GRBs release 
supernova- like energies (~ 10 51 erg). It is thus probable that a major subclass of GRBs is a 
consequence of the collapse of a massive star. 

One of the most attractive scenarios involving massive stars is the collapsar model 
(Woosley 1993; MacFadyen & Woosley 1999). A "collapsar" is a rotating massive star. 
The collapsar model is divided into two classes according to the formation history of the 
black hole. A "type I collapsar" is a failed supernova. The collapse of an iron core leads 
temporarily to a neutron star formation. A black hole, however, is formed quickly within 
a few seconds as a result of the accretion of matter through the stalled shock wave. In the 
mean time, infalling envelope matter is slowed by rotation in the equatorial plane and forms 
a disk. The formation and propagation of relativistic flows from a type I collapsar have been 
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studied numerically in both Newtonian (MacFadyen & Woosley 1999) and special relativistic 
hydrodynamic simulations (Aloy et al. 2000; Zhang, Woosley, & MacFadyen 2003). The 
other model is a "type II collapsar," in which a black hole is formed over a longer period 
of time accompanying a successful supernova. The supernova generates an outgoing shock 
and ejects all the helium and heavy elements outside the neutron star. However, some of the 
post-shock gas fails to reach escape velocity and is pulled back toward the proto-neutron star 
by gravity. When enough gas has fallen back, the neutron star collapses to a black hole. If 
the infalling matter has sufficient angular momentum, it forms a disk at the same time. The 
formation and propagation of relativistic flows from a type II collapsar have been studied 
numerically in Newtonian hydrodynamic simulations (MacFadyen, Woosley, & Heger 2001). 
However, these previous numerical simulations of collapsar models do not fully address the 
outflow mechanism. In these simulations, in fact, the authors estimate the energy of a jet 
assuming neutrino annihilation or magnetohydrodynamic (MHD) process as a source and 
input a jet with that energy from the inner boundary. There have been a lot of studies of 
astrophysical jets with Newtonian (e.g., Uchida & Shibata 1985; Shibata & Uchida 1986, 
1990; Kudoh, Matsumoto, & Shibata 1998; Kato, Kudoh, & Shibata 2002) and relativistic 
MHD simulations (e.g., Koide, Shibata, & Kudoh 1998, 1999). They have fully addressed 
the formation, acceleration, and collimation of jets from accretion disks. We think that the 
formation of relativistic flow from collapsars should also be studied by MHD simulations. 
In fact, several authors (Cameron 2001; Wheeler, Meier, & Wilson 2002) proposed that the 
relativistic jets are generated by the MHD process when the massive star collapses to a 
rapidly rotating neutron star. 

It is suspected that large-scale magnetic fields play an important role in the formation 
of a GRB. Magnetic fields are suitable for extracting energy on the burst time-scale from the 
debris disk in the collapsar scenario. If the field is initially weak, it is likely to be exponentially 
amplified as a result of differential rotation and/or dynamo. Large-scale magnetic fields can 
help guide and collimate the outflow and also contribute to its acceleration (e.g., Usov 1994; 
Thompson 1994; Meszaros & Rees 1997; Katz 1997; Kluzniak & Ruderman 1998; Lyutikov 
& Blackman 2001; Drenkhahn & Spruit 2002; Vlahakis & Konigl 2001, 2003ab). Even if the 
flow is not magnetically driven, the field should be strong enough to account for the observed 
synchrotron emission (Spruit, Daigne & Drenkhahn 2001). If the outflow is largely Poynting 
flux dominated, it can solve the baryon contamination problem. Magnetic fields have been 
favored in this respect for driving GRB outflows. 

The effect of stellar rotation and intrinsic magnetic fields on gravitational collapse of 
massive stars was numerically studied by several authors (LeBlanc & Wilson 1970; Symbal- 
isty 1984; Ardeljan et al. 2000). Symbalisty (1984) computed the magnetorotational core 
collapse of a 15M Q star by numerically solving the 2.5-dimensional MHD equations together 
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with the neutrino transport. The simulations showed the formation of two oppositely di- 
rected, high-density, supersonic jets in the combination of a rapid rotation and a very strong 
dipole magnetic field. Khokhlov et al. (2001) assumed that the jets are generated by a mag- 
netorotational mechanism when a stellar core collapses into a neutron star and simulated 
the process of the jet propagation through the star. Proga et al. (2003) simulated the col- 
lapsar model by using a pseudo- Newtonian MHD code including some essential microphysics 
(physical equation of state [EOS], photodisintegration of nuclei, and neutrino cooling) for 



In this study we perform 2.5-dimensional general relativistic MHD(GRMHD) simula- 
tions of the gravitational collapse of a rotating star with magnetic field as a model for a 
collapsar. The collapsar is in some sense an anisotropic supernova, and it is considered that 
relativistic jets from collapsars are launched by MHD processes in accreting matter and/or 
by neutrino annihilation. We investigate the physics of the formation of jets, the accelera- 
tion force on jets, and the dependence of the acceleration of jets on the initial magnetic field 
strengths and on the initial rotational velocity. We describe the numerical method in our 
simulations in section 2 and present our results in section 3. The summary and discussion 
are given in section 4. 



In order to study the formation of relativistic jets from a collapsar, we use a 2.5- 
dimensional GRMHD code with spherical polar coordinates (r, 0) (i.e., d/d(p = 0, but 
v<p and are nonzero; Koide 2003. The method is based on the general relativistic formula- 
tion of the conservation laws of particle-number and energy-momentum, Maxwell equations, 
and Ohm's law with no electrical resistance (ideal MHD condition) on a curved space-time 
(Weinberg 1972; Koide et al. 1998, 1999; Koide 2003). 

The space-time (x°, x 1 , x 2 , x 3 ) = (ct, x 1 , x 2 , x 3 ) is described by the metric g^, where 
the line element ds is given by (ds) 2 = g ^dx^ dx u . Here c is the speed of light. The basic 
equations of GRMHD in four-dimensional space-time are 
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Numerical Method 



2.1. Basic Equations 
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d„F vX + d u F x ^ + d\F^ u = 0, (3) 

V,F^ = -pt r, (4) 

where V u is the covariant derivative, \\g\\ is the determinant of g^ v (as a matrix), and 
= \g x ° + f|f + ^) are the Christoffel symbols. U v and J" = (cp e , J 1 , J 2 , J 3 ) 

are four- velocity and four-current density, respectively (p e is the electric charge density). 
The general relativistic energy momentum tensor T^ u is given as 

7*" = p<T + (e itrt + p)UW + F%F VC - ^F Xk F Xk (5) 

where _F Miy is the electromagnetic field-strength tensor, F^ v = d^A v — d^A^ and = 
($ e /c, A 1 , A 2 , A 3 ) is the four-vector potential (<£> e is the electro-static potential). The elec- 
tric field Ei and the magnetic field Bi are given by Ei = cF i0 (i = 1,2,3) and B x = F 23 , 
B2 = -F31, B 3 = F12, respectively. The electric field Ei and the magnetic field Bi are normal- 
ized as Ei = E*/y/JI and Bi = B* / \[Jl, respectively, where \x is the magnetic permeability 
and quantities with asterisks are given in MKSA units (SI units). The scalars p, p, and t int 
are proper mass density, proper pressure, and proper internal energy density, respectively. 
We assume mainly for numerical simplicity that matter can be described as an ideal gas; 
hence, = pc 2 +p/(T — 1), where T is the specific heat ratio although we know that the 
gamma law EOS is not very appropriate for the gravitational collapse of massive stars. We 
further assume that T = 5/3 in the simulations for numerical reasons although it is more 
appropriate to adopt T = 4/3 for the gas of current interest. In addition to the equations, 
we assume the infinite electric conductivity condition 

F^U U = 0. (6) 

Then, equations (l)-(3) together with equation (6) are a closed system of equations. Equation 
(4) is used only to calculate the current density J M . 

If we assume that the off-diagonal elements of the metric g^ vanish, 

g»u = (/i ^ v), (7) 

and we use the notation 

#00 = -hi, g n = hi, g 22 = h\, g 33 = h\, (8) 
then the line element can be written as 

3 

(dsf = g^dafdx" = -h 2 (cdt) 2 + J2 h\(dx 1 ) 2 . (9) 

i=i 
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We employ the 3+1 formalism (Thorne, Price, & Macdonald 1986; Koide et al. 1999; 
Koide 2003). The velocity of light is written explicitly. In the non-relativistic limit (c — > oo), 
they reduce to the standard Newtonian MHD equations. 

In the GRMHD code, a simplified total variation diminishing (TVD) method is employed 
(Davis 1984). This method is useful because it requires only the maximum speed of physical 
waves but not the eigenvalues of the Jacobian of the equations. 

De Villiers & Hawley (2003) and Gammie, McKinney, & Toth (2003) have developed 
a new GRMHD code independently. De Villiers & Hawley (2003) used a non- conservative 
scheme that is based on techniques developed by Hawley, Smarr, & Wilson (1984) for ax- 
isymmetric hydrodynamics around black holes. On the other hand, Gammie, McKinney, & 
Toth (2003) used a conservative, shock-capturing scheme. They have performed some one- 
and two- dimensional tests and have obtained consistent results with each other. 

We do not consider the evolution of the metric because the accreted mass is sufficiently 
small in the time-scale of the simulations. The current calculation focuses on GRMHD 
with simplified microphysics. Although photodisintegration of bound nuclei and cooling by 
neutrino emission (thermal and especially Urea processes) have been shown to be of critical 
importance for the dynamics of collapsar accretion disks (Popham, Woosley & Fryer 1998; 
MacFadyen & Woosley 1999), we ignore these effects in the current calculations for simplicity. 
It is hoped that future GRMHD simulations will include essential microphysics (neutrino 
cooling, physical EOS, photodisintegration) to realistically simulate collapsar accretion disks. 

2.2. Initial Condition 

As for the initial model, we have the collapsar model in mind. In principle, we should 
start calculations from a realistic progenitor model with rotation and magnetic field. Al- 
though some progress has been made lately on the precollapse evolution of massive stars 
with rotation and magnetic field (Heger et al. 2003), their models have a lot of uncertainty 
yet. Hence, we take the following pragmatic approach: Core collapse is assumed to fail in 
this star, and a few M & black hole is inserted by hand into the calculation. The simulations 
presented in the paper follow the accretion of gas into a black hole assumed to have formed 
before the calculation begins. Furthermore, We use Bruenn's realistic one-dimensional su- 
pernova model (Bruenn 1992) as a rough guide to what the density structure exterior to the 
black hole might be. The gas in Bruenn's model must have collapsed before the present 
calculation began to form the black hole that is assumed to be present. However, we use 
the scaling in that model as a guide for the scale-free stellar structure. We employ only the 
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profiles of the density, pressure, and radial velocity as our initial condition. In this way, we 
can discuss generic features of the dynamics. 

We consider a non-rotating black hole. The metric is given by 

900 = 1 9ll= i 1 ~~i) ' 922 = R2 > 933 = R 2 sme, (10) 

where r$ = 2GM Q /c 2 is the Schwarzschild radius. We add the effect of stellar rotation and 
intrinsic magnetic field in the simulations as the initial condition. The initial rotational 
velocity distribution is assumed to be a function of the distance from the rotation axis, 

r = Rsm9, (11) 

only: 

2 

v * = Vo ^hi r - (12) 

Here, v is a model parameter for rotational velocity. We fix x = 100 r s in this paper. The 
initial magnetic field is assumed to be uniform and parallel to the rotational axis. This is 
known as the Wald solution for a non-rotating black hole: Br = Bq cos 9, Be = —aBo sin 9. 
Bq is the overall magnetic field strength and a = goo is the lapse function. Figure 1 shows the 
schematic picture of our simulation. The distributions of various physical quantities (density 
p, pressure p, each velocity component v r ,v^,v z , each magnetic field component B r , B^, B z , 
and plasma beta f3 = Pgas/Vmag) on the equatorial plane for the initial state are shown in 
Figure 2. In the distribution of density and pressure, we can see that the weak shock is 
sitting at about 18 rs- We input the stellar rotation up to 18 r s because the accreted mass 
for stellar matter in the pre-shock region is sufficiently small in the time scale. We neglect 
rotation in the pre-shock region for numerical simplicity. The rotation is almost uniform 
(see Eq. [12]). This initial rotation profile is similar to the previous simulations of rotational 
core collapse (Monchmeyer & Miiller 1989; Yamada & Sato 1994). 

We want to emphasize that our simulations are scale free. The normalization units are 
summarized in Table 1. We use typical values for normalizations. If we set a black hole mass 
Mbh ~ 3M , the Schwarzschild radius rs is about 8x 10 5 cm and the time unit rg(= rs/c) is 
3.0 x 10~ 5 s. The units of magnetic field strength and pressure depend on the normalization 
of the density. If we take, for example, the density unit po — 10 10 g/cm 3 , the magnetic 
field strength unit is 3.0 x 10 14 G and the pressure unit is P = 10 31 dyn/cm 2 . The models 
computed in this paper are summarized in Table 2. These units are the values at the last 
stable Keplerian orbit of a non-rotating black hole on the equatorial plane \{r,z) = (r ,0)]. 
The variables with subscript "0" represent the values at the last stable orbit throughout this 
paper. 
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We also use the following nondimensional parameters: 

KO 



E t h - Trf-, (13) 



V 2 

E =-M (U) 

V K0 



E ro t = JJ2~, (15) 



V 2 

V K0 



where V s0 is the relativistic sound speed, 



p c 2 + r> /(r - 1) 



Vao is the Alfven velocity, 



V A0 = c — B ° (17) 

v/p c 2 + r> /(r - 1) + b 2 ' 

and Vko is Keplerian velocity, 

V K0 = c/[2(r /r s -l)} 1 / 2 . (18) 

They are all evaluated at the last stable orbit. The values of the parameters adopted in this 
paper are summarized in Table 2. 

The simulations are done in the region 2r s < R < 60r 5 , < 9 < n/2, with 120 x 120 
mesh points. We assume axisymmetry with respect to z-axis and mirror symmetry with 
respect to the equatorial plane. We employ a free boundary condition, which waves, fluids 
and magnetic fields can pass through freely, at R = 2r s and R = 60r s (see Koide et al. 
1999). 



3. Numerical Results 

3.1. Formation of Jet 

We shall first discuss the results of the case A3 (Table 2; from now on referred to as 
the standard case). The model parameters are B = 0.05, v = 0.01, E th = 2.36 x 10~ 3 , 
Emag = 1-68 x 10~ 3 , and E rot = 5.36 x 10~ 2 . Figure 3 shows the time evolutions of density. 
The stellar matter falls onto the central black hole at first. This collapse is anisotropic as a 
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result of the effects of rotation and magnetic field. The accreting matter falls more slowly on 
the equatorial plane than on the rotational axis. The matter piles up on the equatorial plane, 
and a disklike structure is formed near the central black hole. At t/rs = 120, the disk density 
is about 2 orders of magnitude higher than the density of matter above the disk. Since the 
magnetic field is frozen into the plasma, it is dragged by the accreting matter and deformed 
to be more radial. At about t/rs — 60, a shock wave is formed near the central black hole and 
propagates outward anisotropically, because of the geometry of the magnetic field. A vortex 
pattern is generated inside the shock front, and moves outwardly with the propagating shock 
front, and produces a jetlike outflow. As this jetlike outflow moves up along the magnetic 
field lines, it becomes collimated. The simulation is terminated at t/r s = 175, at which time 
the jet is still propagating outward. 

The jetlike outflow is generated behind the shock wave. Figure 4 shows the time evo- 
lutions of plasma beta (= P gas /Pmag) and toroidal magnetic field for the standard case A3. 
The shock wave expands with twisted magnetic fields, and the outflow has high magnetic 
pressures. The jetlike outflow formed in the simulation is thus magnetically driven. 

The magnetic field plays an important role also in the collimation. The plasma beta is 
initially low in the whole region, indicating that the magnetic pressure is dominant over the 
gas pressure (see Fig. (4)). Not only the pinching force of the toroidal magnetic field but also 
the geometry of the poloidal magnetic field (i.e. poloidal magnetic pressure) plays a crucial 
role. 



3.2. Mechanism of Shock 

We now investigate in more detail the mechanism of shock formation in the simulation. 
Figure 5 shows the time evolution of the spatial distributions of plasma beta (f3 = Pgas/Pmag) 
and toroidal magnetic field near the black hole in the standard case A2. At t/rs — 60.0, the 
accreting stellar matter rotates faster than earlier around the central black hole. As time 
goes on, the magnetic field is more and more twisted as a result of the differential rotation of 
accreting matter after the poloidal magnetic field is deformed to be more radial by the radial 
accretion. The toroidal magnetic field is amplified significantly near the central black hole. 
The amplified magnetic field expands outward as Alfven waves and launches an outgoing 
shock wave. We can see that the expansion of the toroidal magnetic field takes place at the 
same position as the shock front. The plasma beta {(3 = Pgas/Pmag) inside the shock wave 
is low {(5 < 1), implying that the shock wave is generated mainly by the magnetic field. 

Figure 6 and 7 show the time evolutions of various physical quantities (density, gas 
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and magnetic pressures, magnetic field, velocity, and forces) on the cylindrical surface with 
x/r s = 4. It shows the propagation of the outgoing shock. It is also found that when 
the shock wave is generated near the central black hole, the magnetic field is amplified by 
about a factor of 6 from the initial strength of the magnetic field near the black hole. The 
magnetic force F mag is comparable to the centrifugal force F cen . These are about an order 
of magnitude higher than the pressure gradient force F p . This confirms the claim that the 
magnetic force, together with the magnetocentrifugal force produces the shock wave rather 
than the pressure gradient force. 



3.3. Properties of Jet 

We discuss here the properties of the jet found in our simulation. Figure 8 shows the 
distributions of various physical quantities along the jetlike outflow at t/rs = 175. By this 
time, the shock wave has already passed through the upper boundary [zjr s = 30) of this 
figure. 

It is easily seen in the velocity distribution that the jet has a mildly relativistic velocity, 
~ 0.3 c (the poloidal velocity is ~ 0.1 c). It is larger than the sound velocity, as well as 
the Keplerian velocity. It also clearly exceeds the escape velocity at that point, and the jet 
is likely to get out of the stellar remnant. The toroidal velocity is larger than the poloidal 
velocity in the jet-forming region. The density of the jet is about twice as high as the density 
of the surrounding matter, ~ 3 x 10 8 g/cm 3 if the density unit is po — 10 10 g/cm 3 . The 
toroidal magnetic field is the dominant component of the magnetic field in the jet-forming 
region. The total magnetic field strength there is ~ 3.0 x 10 13 G. 

On the other hand, we see in the energy distribution that the magnetic energy is compa- 
rable with the rotational energy and larger than the gravitational energy, which is consistent 
with the previous statement that the jet is accelerated by the magnetic and centrifugal forces. 

Figure 9 shows the time variation of the mass flux of the jet and the accretion rate, 
the kinetic energy and the Poynting flux, and the jet and Alfven wave components of the 
Poynting flux at z/r s ^ 12 for the standard case A2. Here, the jet and Alfven components 
are defined as 

B 2 

SEMjet = -7—V z , (19) 
47T 

SEMAlf = —J—^ v <t>- (20) 

The Poynting flux is an order of magnitude larger than the kinetic energy flux. In the 
Poynting flux, the Alfven wave component dominates over the jet components. This means 
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that the Alfven wave transports more energy outward. Using the simulation data, we can 
estimate that the kinetic energy of the jet, Ej et , is ~ 10 49 ergs if the density of the jet 
is 10 10 g/cm 3 . This is 2 orders of magnitude lower than the standard energy of GRBs (~ 
10 51 ergs). If Poynting flux is considered, the total energy becomes ~ 10 50 ergs, still less 
than the energy of GRBs. 



3.4. Dependence on the Initial Magnetic Field Strength 

The dependence of the jet properties on the initial magnetic field strength has also been 
investigated. Figure 10 shows the density and plasma beta distributions for case A3 in which 
the initial magnetic field is greater than the standard case A2 (B = 0.07, Vq = 0.01). The 
jetlike outflow is weaker and fainter than in the standard case A2. The maximum poloidal 
velocity of the jet is less than 0.05 c, although the jet velocity is still larger than the escape 
velocity. Thus, even in this case, a jet will emerge. The shock wave itself, as well as the 
magnetic twist, is also weaker. 

The opposite is true for weaker magnetic fields. Shown in Figure 11 is model Al (Bo = 
0.03, vq = 0.01). The jetlike outflow is clearly stronger than in the standard case A2. The 
maximum poloidal velocity of the jet becomes about 0.1c. The total velocity exceeds the 
escape velocity inside the shock. The magnetic twist is significantly stronger, and as a result 
the propagating shock wave is stronger as well. 

Figure 12 summarizes these dependenices on the initial magnetic field strength. As the 
initial magnetic field strength increases, the jet velocity (v z , v^) increases and the magnetic 
twist decreases. However, for stronger magnetic field (B > 0.055) the jet velocity decreases 
with increasing B and the magnetic twist still continues to decrease. In order to produce a 
strong shock, the magnetic field has to be twisted significantly so that it can store enough 
magnetic energy. If the initial magnetic field is strong, the magnetic field cannot be twisted 
significantly because Alfven waves propagate as soon as the magnetic field is twisted a little 
bit. As a result, the jet velocity does not rise up and the magnetic twist remains weak. 
Therefore, weaker initial magnetic fields are favorable for a stronger jet. 

Let us now consider the physical reason of the above results. In the Newtonian case, 
the time evolution of the toroidal magnetic field is given as 



dt 

where 



0D " uB p , (21) 



1 / X 

to ~ — (22) 

trot 
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and t rot is the rotation time scale. From this we obtain 



n wt. (23) 

-Dp 



In this situation, the comparing time scale is the propagation time scale of the Alfven wave 



L L^Fp 

r ' ^ 

VAz Hz 

where L is the typical length and va z is the vertical component of the Alfven wave velocity 
vector. Substituting equation (24) for equation (23), we have 

B, ^ uWWp k J_ (25) 
Bp B z B z 



In this simulation, B z ~ B . Therefore, 

B d 1 



f-NJ 



-Bp -Bq 



(26) 



This can explain the dependence of the magnetic twist on the initial magnetic field strength, 
especially in the case of the Alfven wave component (see Fig. 11c, dashed line). 

The upward motion of the fluid behind the twist wave front is induced by the J x B 
force. If we neglect other forces, the equation of motion for the fluid element in the z-direction 
becomes 

P^-vfUV (27) 



dt \8n, 
which can be rewritten as 

«. --*-(?)■ P8) 

p z \8n J 

In this situation the time scale is determined by the propagation time scale of the Alfven 
wave. Thus, t/z ~ 1/vaz- Since in the region near the wave front, B p ~ B z ~ B , equation 
(28) can be rewritten as 

1 1 / B 2 \ 

' ! / X ocB \ (29) 



v z ~ 

pv Az \8n 

This explains the dependence of v z on B for B > 0.055 in Figure 11a; i.e., as the ini- 
tial magnetic field strength increases, the vertical component of the jet velocity decreases. 
However, the jet velocity has a limit when the initial magnetic field strength becomes even 
weaker. It can be understood as follows: the energy transport by the magnetic field is the 
Poynting flux 

c B z Bq 

— E x B ~ -Zv g + —t^vq. 30 

47T 47T 471 V 
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The first term on the right-hand side of equation (30) is the jet component, whereas the 
second term is the Alfven wave component. This equation also can explain the dependence 
of the jet velocity on the initial magnetic field strength. When the initial magnetic field 
strength becomes stronger, B z increases. This means that the Alfven wave component takes 
away more energy than the jet component. Thus, a less magnetic energy is converted to 
the kinetic energy of the jet, so that the jet becomes weaker. When the initial magnetic 
field strength becomes weaker, the magnetic twist becomes stronger and increases. Thus, 
much of the magnetic energy is converted to the kinetic energy of the jet. However, the 
converted magnetic energy has a limit. When the magnetic twist (B p /B^) of the Alfven 
wave becomes 1, it becomes a limit. This threshold is determined by the balance between 
the jet component and the Alfven wave component 

B\ B z B,p 

4^ = 17*' (31) 

which can be rewritten as 

^ = V - + ~ 1. (32) 

B z v z y ' 

Thus, we can understand that when the magnetic twist in the Alfven wave becomes 1, the 
jet velocity becomes maximum. 

On the other hand, when the initial magnetic field is zero (case D), the shock wave is 
not generated and the jetlike outflow is not formed either. This clearly demonstrates that 
the generation of the shock wave, as well as the jetlike outflow, requires magnetic fields. 



3.5. Dependence on the Initial Rotational Velocity 

Next we show the dependence of the jet properties on the initial rotational velocity. 
Figure 13 shows the results (density and plasma beta distributions) of the fast rotation case 
B2 (Bo = 0.05, vq = 0.015). The jetlike outflow is stronger and wider than in the standard 
case A2. The maximum poloidal velocity of the jet is about 0.1 c. The jet velocity exceeds 
the local escape velocity in the whole region. The magnetic twist is stronger accordingly. 

The slow rotation case Bl (Bo = 0.05, vq = 0.005) is shown in Figure 14. The jetlike 
outflow is weaker and fainter. The maximum poloidal velocity of the jet is less than 0.01 
c. The jet velocity is barely larger than the escape velocity in some regions. The magnetic 
twist is weaker correspondingly. 

Figure 15 shows the dependence of various physical quantities of the jet on the initial 
rotational velocity. As the initial rotational velocity becomes faster, the jet velocity (v z , v^) 
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becomes faster and the magnetic twist becomes stronger up to a certain value (t> = 0.01). 
For faster initial rotational velocity (v > 0.01), the vertical component of the jet velocity (v z ) 
and the magnetic twist are almost constant and the toroidal component of the jet velocity 
(v,/,) becomes a little bit slower. If the initial rotation is sufficiently fast, the magnetic field 
is twisted significantly and stores enough energy to produce a fast jet. However, the stored 
magnetic energy has a limit. It is determined by the competition between the propagation 
time of Alfven waves and the rotation time of the disk (or the twisting time). This can be 
understood by the same reason as the dependence on the initial magnetic field strength. The 
limit is determined when the magnetic twist of the Alfven wave becomes 1. Provided with 
the results of this and the previous sections, the magnetic field strength is the dominant 
factor to determine the maximum jet velocity. It should be noted, however, that when the 
initial rotational velocity is zero (case C), the shock wave is not generated and the jetlike 
outflow is not formed either. Hence, the magnetic twist is crucial to produce the shock wave 
and the jetlike outflow. 



4. Summary and Discussion 

We have studied the generation of a jet from gravitational collapse of a rotating star 
with magnetic fields by using the 2.5-dimensional GRMHD simulation code. Our results are 
summarized as follows: 

1. When the stellar matter falls onto the black hole, the collapse is anisotropic because 
of the effect of rotation and magnetic field and a disklike structure is formed near the 
black hole. A shock wave is launched near the black hole and propagates outward. 
The jetlike outflow is produced inside the shock. The maximum velocity of the jet is 
mildly relativistic (~ 0.3 c). This result is consistent with the pseudo-Newtonian case 
(Proga et al. 2003). The kinetic energy of the jet is ~ 10 49 ergs. 

2. The shock wave is generated mainly by magnetic and centrifugal forces. The frozen 
magnetic field is twisted and amplified by the rotation of stellar matter particularly 
near the black hole. It starts to expand outwardly with Alfven waves and leads to the 
production of the shock wave. 

3. The acceleration of the jet is also driven by the magnetic pressure and centrifugal 
forces. The jet is collimated in the course of propagation by the pinching force of the 
toroidal magnetic field and the geometry of the poloidal magnetic field. 

4. The magnetic twist is large in the jet-forming region. The jet has a helical structure. 
This is similar to other astrophysical jets formed by MHD processes. 
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5. The Poynting flux is an order of magnitude larger than the kinetic energy flux. Alfven 
waves as the Poynting flux transports more energy outward than the jet. 

6. As the initial magnetic field strength increases, the jet velocity increases and the mag- 
netic twist decreases. However, for stronger magnetic field the jet velocity decreases 
with increasing initial magnetic field strength and the magnetic twist goes on decreas- 
ing. The dependence of the jet properties on the initial rotational velocity is similar 
to the dependence on the initial magnetic field strength. As the initial rotational ve- 
locity becomes faster, the jet velocity becomes faster and the magnetic twist becomes 
stronger up to a certain value. For faster initial rotational velocity, the jet velocity 
and the magnetic twist are almost constant. The magnetic energy is converted to the 
kinetic energy of the jet. However, the converted magnetic energy has a limit. It is 
determined by the competition between the propagation time of Alfven waves and the 
rotation time of the disk (or the twisting time). Hence the jet velocity is saturated at 
some point. 

The maximum jet velocity in the current simulations was about 0.3 c and the typical 
kinetic energy of the jet was ~ 10 49 ergs . The jet is too slow and too weak for the jet 
of GRBs. However, the baryon mass in the jet is small. This gives good results for the 
GRB model. We have to consider other acceleration mechanisms. We expect that disk 
jets formed from an accretion disk become faster than the jet in the current simulations 
because they can release more gravitational energy and have enough energy to be applied 
to the GRB model. We believe that, if the current simulations can be continued for longer 
times, a disk-jet may be formed by the MHD processes. If the jet contains enough energy 
to convert the kinetic energy of the jet, then the break-out of the jet through the stellar 
surface is a possibility. When the jet goes through the stellar surface, the strong density 
gradient may accelerate the jet. In fact, some authors (Aloy et al. 2000; Zhang, Woosley, & 
MacFadyen 2003) have shown numerically that a significant acceleration of the jet occurs and 
the terminal Lorentz factor becomes as high as T ~ 50. Since they used special relativistic 
hydrodynamics codes and neglected magnetic fields altogether, we think that it is important 
to simulate the propagation of the jet outside the stellar surface by a GRMHD code properly 
evaluating the importance of magnetic fields for the dynamics and further propagation of 
the jet. 

On the other hand, our results may apply to baryon-rich outflows associated with the 
so-called failed GRBs. It is supposed to be a high baryon-load fireball attaining mildly 
relativistic velocities and not producing GRBs. Such failed GRBs may have event rates higher 
than successful GRBs (Woosley et al. 2002; Huang, Dai, & Lu 2002). Some failed GRBs can 
be observed as "hypernovae." SN 2002ap may be such an event. It shows broad-line spectral 
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features like the famous SN 1998bw (Kinugasa et al. 2002), yet it is not associated with 
a GRB. Kawabata et al. (2002) proposed a jet model for SN 2002ap consistent with their 
spectropolarimetric observations. Totani (2003) estimated that the jet should be moving at 
about 0.23 c in the direction perpendicular to the line of sight and that the kinetic energy 
of the jet was ~ 5 x 10 50 ergs . 

We have assumed a uniform global magnetic field in the simulations. However, this 
assumption may be unrealistic. From observations, it is inferred that rotating compact stars 
in general have a dipole-like magnetic field. Hence, it may be more likely that the rotating 
stars collapse with a dipole-like magnetic field. On the other hand, it is possible that the 
magnetic field near the central black hole is radial because of radial accretion of matter. We 
will study these possibilities in a forthcoming paper. There are a lot of works on astrophysical 
jets using non-relativistic MHD simulations with a dipole magnetic field configuration (e.g., 
Hayashi, Shibata, & Matumoto 1996; Goodson, Winglee, & Bohm 1997; Goodson, Bohm, 
& Winglee 199; Goodson & Winglee 1999; Romanova et al. 2002; Kato et al. 2003). These 
simulations show a high- velocity wind and lead the ejection of hot plasmoids. If the plasmoids 
are ejected intermittently, they will collide with each other, forming "internal shocks," which 
may explain the time variability of GRBs (Shibata & Aoki 2003; Aoki, Yashiro, & Shibata. 
2003). 

Y. M. appreciates many helpful conversations on GRBs with S. Aoki, M. Uemura, 
and T. Totani. He also thanks K. Watarai, Y. Kato, K. Uehara, K. Nishikawa, and T. 
Haugboulle for useful discussions. This work was partially supported by Japan Science and 
Technology Cooperation (ACT-JST), Grant-in- Aid for the 21st Century COE "Center for 
Diversity and Universality in Physics" and Grants-in-Aid for the scientific research from the 
Ministry of Education, Science, Sports, Technology, and Culture of Japan through 14079202, 
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Observatory, Japan (yyml7b), and partly on the Alpha Server ES40 at the Yukawa Institute 
for Theoretical Physics of Kyoto University, Japan. 

REFERENCES 

Aloy, M. A., Muller, E., Ibanez, J. M., Marti, J. M., & MacFadyen, A. I. 2000, ApJ, 531, 
L119 

Aoki, S., Yashiro, S, & Shibata, K. 2003, preprint 

Ardeljan, N. V., Bisnovatyi-Kogan, G. S., & Moiseenko, S. G. 2000, Ap&SS, 274, 389 



-17- 

Bloom, J. S., Kulkarni, S. R., & Djorgovski G. 2002, A J, 123, 1111 
Bloom, J. S., et al. 2002, ApJ, 572, L45 

Bloom, J. S., Frail, D. A., & Kulkarni, S. R. 2003, ApJ, 594, 674 

Bruenn, S. W. 1992, in Neuclear Physics in the Universe, ed. M. W. Guidry & M. R. Strayer 
(Philadelphia: Institute of Physics Pub.), 31 

Cameron, A. G. W. 2001, ApJ, 562, 456 

Davis, S. F. 1984, NASA Contractor Rep. 172373, ICASE Rep., No. 84-20 
De Villiers, J. -P. & Hawley, J. F. 2003, ApJ, 589, 458 
Drenkhahn, G. & Spruit, H. C. 2002, A&A, 391, 1141 
Fishman, G. J. & Meegan, C. A. 1995, ARA&A, 33, 415 
Frail, D., et al. 2001, ApJ, 562, L55 

Gammie, C. F., McKinney, J. C, & Toth, G. 2003, ApJ, 589, 444 
Garnavich, P. M., et al. 2003, ApJ, 582, 924 

Goodson, A. P., Winglee, R. M., & B6hm, K. H. 1997, ApJ, 489, 199 

Goodson, A. P., B6hm, K. H., & Winglee, R. M. 1999, ApJ, 524, 142 

Goodson, A. P. & Winglee, R. M. 1999, ApJ, 524, 159 

Hawley, J. F., Smarr, L. L., & Wilson, J. R. 1984, ApJS, 55, 211 

Hayashi, M. R., Shibata, K., & Matumoto, R. 1996, ApJ, 468, L37 

Heger, A., Woosley, S. E., Langer, N, & Spruit, H. C. 2003, preprint (astro-ph/0301374) 

Hjorth, J., et al. 2003, Nature, 423, 847 

Huang, Y. F., Dai, Z. C, & Lu, T. 2002, MNRAS, 332, 735 

Iwamoto, K., et al. 1998, Nature, 395, 672 

Kato, S. X., Kudoh, T., & Shibata, K. 2002, ApJ, 565, 1035 

Kato. Y., Mineshige, S., & Shibata, K. 2003, ApJ, accepted (astro-ph/0307306) 



Katz, J. I. 1997, ApJ, 490, 633 

Kawabata, K. S., et al. 2002, ApJ, 580, L39 

Kinugasa, K., et al. 2002, ApJ, 577, L97 

Khokhlov, A. M., Hoflich, P. A., Oren, E. S., Wheeler, J. C, Wang, L., & Chtchelkanova, 
A. Yu. 1999, ApJ, 524, L107 

Klebesadel, R. W., Strong, I. B., & Olson, R. A. ApJ, 182, L85 

Kluzniak, W. & Ruderman, M. 1998, ApJ, 505, LI 13 

Koide, S., Shibata, K., & Kudoh, T. 1998, ApJ, 495, L63 

-. 1999, ApJ, 522, 727 

Koide, S. 2003, Phys. Rev. D, 67, 104010 

Kudoh, T., Matsumoto, R., & Shibata, K. 1998, ApJ, 508, 186 
LeBlanc, J. M. & Wilson, J. R. 1970, ApJ, 161, 541 
Lyutikov, M. & Blackman, E. G. 2001, MNRAS, 321, 177 
MacFadyen, A. I. & Woosley, S. E. 1999, ApJ, 524, 262 
MacFadyen, A. I., Woosley, S. E. & Heger, A. 2001, ApJ, 550, 410 
Meszaros, P. & Rees, M. J. 1997, ApJ, 482, L29 
Meszaros, P. 2002 ARA&A, 40, 137 

Munchmeyer, R. M. & Miiller, E. 1989, in NATO ASI Series C262, Timing Neutron Stars, 
ed. H. Ogelman & E.P.J. Van den Heuvel (New York: ASI), 549 

Panaitescu, A. & Kumar, P. 2001, ApJ, 560, L49 

Popham, R., Woosley, S. E., & Fryer, C. 1999 ApJ, 518, 356 

Piran, T. 1999, Phys. Rep., 314, 575 

Price, P. A., et al. 2003, Nature, 423, 844 

Proga, D., MacFadyen, A. I., Armitage, P. J., & Begelman, M. C. 2003, ApJ, 599, L5 



-19- 



Romanova, M. M., Ustyugoya, G. V., Koldoba, A. V., & Lovelace, R. V. E. 2002, ApJ, 578, 
420 

Sari, R., Piran, T., & Halpern, J. P. 1999, ApJ, 517, L17 
Shibata, K. & Aoki, S. 2003, preprint (astro-ph/0303253) 
Shibata, K. & Uchida, Y. 1986, PASJ, 38, 631 

-. 1990 ,PASJ, 42, 39 
Spruit, H. C, Daigne, F., & Drenkhahn, G. 2001, A&A, 369, 694 
Stanek, K. Z., et al. 2003, ApJ, 591, L17 
Symbalisty, E. M. D. 1984, ApJ, 285, 729 
Thompson, C. 1994, MNRAS, 270, 480 

Thorne, Kip S., Price, R. H., & Macdonald, D. A. 1986, Membrane Paradigm (Yale Univer- 
sity Press, New Haven and London) 

Totani, T. 2003, ApJ, in press (astro-ph/0303621) 

Uchida, Y. & Shibata, K. 1985, PASJ, 37, 515 

Uemura, M., et al. 2003, Nature, 423, 843 

Usov, V. V. 1994, MNRAS, 267, 1035 

van Paradijs, J., Koubelitou, C, & Wijers, Ralph A. M. J. 2000, ARA&A, 38, 379 
Vlahakis, N. & Konigl, A. 2001, ApJ, 563, L129 

-. 2003a, ApJ, 596, 1080 

-. 2003b, ApJ, 596, 1104 
Weinberg, S. 1972, Gravitation and Cosmology (wiely, New York) 
Wheeler, J. G, Meier, D. L., & Wilson, J. R. 2002, ApJ, 568, 807 
Woosley, S. E. 1993, ApJ, 405, 273 

Woosley, S. E., Eastman, R. G., & Schmidt, B. 1999, ApJ, 516, 788 
Woosley, S. E., Zhang, W. & Heger, A. 2002, preprint(astro-ph/0211063) 



-20- 

Yamada, S. & Sato, K. 1994, ApJ, 434, 268 

Zhang, W., Woosley, S. E., & MacFadyen, A. I. 2003, ApJ, 586, 356 



This preprint was prepared with the AAS IATgX macros v5.2. 



-21- 



z a 

magnetic licit! lin« simulation box 



ir^L 1 boundary 




L^LLjLuiriaJ mirror symmetry 



Fig. 1. — Schematic picture of our simulations. Dashed lines show the simulation box. The 
gray region shows post-shock matter. The black region shows the black hole. 
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Fig. 2. — Various physical quantities on the equatorial plane z = at the initial time in 
the standard case A2 (B = 0.05, vq = 0.01). (a) Density p (solid line) and gas pressure P 
(dashed line), (b) Components of velocity v r (solid line), (dashed line), and v z (dot-dashed 
line); vk (double- dot- dashed line) is the Keplerian velocity, (c) Components of magnetic field 
B r (solid line), (dashed line), and B z (dot-dashed line), (d) Plasma beta (3 = P gas / 1 P ma g- 
The stellar matter falls toward the black hole and rotates like a rigid body around the black 
hole. 
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Fig. 3. — Time evolution of the density for case A2 (B = 0.05, vq = 0.01). The color scale 
shows the value of the logarithm of density. The white curves depict magnetic field lines. 
Arrows represent the poloidal velocities normalized by the light velocity. The central black 
region corresponds to the inner boundary at r/r$ = 2. The distance and the time are given 
in units of r s and r s = r s /c, respectively, (a) Initial condition. The stellar matter falls and 
rotates around the black hole, (b) Condition at t/r s = 60. Accreting matter forms a disklike 
structure, (c) Condition at tjr s = 110. The shock wave is generated near the black hole. 
(d) Condition at t/r s = 175. The shock wave propagates outward. The jetlike outflow is 
produced inside the shock wave. 
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Fig. 4. — The time evolution of the plasma beta (j3 = P gas /Pmag) distribution for the 
standard case A2. The color scale shows the value of the logarithm of plasma beta. The white 
contour plots are the toroidal component of the magnetic field. Arrows show the poloidal 
velocities normalized by the light velocity. The central black region corresponds to the inner 
boundary. The distance and the time are given in units of r s and r s = r s /c, respectively. 
The amplified magnetic field expands outward and collimates the jetlike outflow. 
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Fig. 5. — Time evolution of the plasma beta {(3 = P ga s/ Pmag) distribution near the central 
black hole for case A2. The color scale indicates the value of the logarithm of plasma beta. 
The white contour plot shows the toroidal component of the magnetic field. Arrows present 
the poloidal velocities normalized by the light velocity. The distance and the time are given 
in units of rs and t$ = rs/c, respectively, (a) Condition at t/rs = 45. The toroidal 
magnetic field is amplified near the black hole. The plasma beta is low in the whole region. 
(b) Condition at t/r s = 60. The shock wave is generated near the black hole, (c) Condition 
at t/r s = 80. The shock wave propagates outwards with amplified toroidal magnetic fields. 
(d) Condition at t/r s = 110. The jetlike outflow is generated inside the shock wave. The 
plasma beta is low in the post-shock region. 
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Fig. 6. — Time evolutions of the various physical quantities on the cylindrical surface of 
x = 4.0r s . (a, b, c) Density p (solid line), gas pressure P gas (dashed line), and magnetic 
pressure P mag (dot-dashed line), (d, e, f) Gravitational force F gr (solid line), gas pressure 
gradient force F p (dashed line), magnetic force F mag (dot-dashed line), and centrifugal force 
F cen (double-dot-dashed line). The times are the same as in Figure 5. These panels show the 
propagation of the shock (see arrows). In the post-shock region, density, gas pressure, and 
magnetic pressure increase. In (f), the magnetic and centrifugal forces are larger than the 
gravitational force. 
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Fig. 7. — Time evolution of the various physical quantities on the cylindrical surface of 
x = 4.0r,5. ( a , b, c) show the components of velocity, v r (solid line), (dashed line), and 
v z (dot-dashed line), (d, e, f) Components of magnetic field, B r (solid line), (dashed 
line), and B z (dot-dashed line). The times are the same as in Figure 5 and 6. The amplified 
magnetic field expands outward. Accreting matter is pushed back and the jetlike outflow is 
generated. 
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Fig. 8. — Various physical quantities on the surface of x/rs = 12 at t/rs = 175 for case A2. 
(a) Density p (solid line), gas pressure P gas (dashed line), and magnetic pressure P mag (dot- 
dashed line), (b) Components of velocity, v r (solid line), (dashed line), and v z (dot-dashed 
line), (c) Matter velocity v tot (solid line), Alfven velocity va (dashed line), sound velocity 
cs (dot-dashed line), and Keplerian velocity vk (double- dot- dashed line), (d) Components 
of magnetic field, B r (solid line), (dashed line), and B z (dot-dashed line), (e) Ratio 
of toroidal to poloidal magnetic field components B^/B p . (f) Gravitational energy E gr 
(solid line), thermal energy E th (dashed line), magnetic energy E mag (dot-dashed line), and 
rotational energy E rot (double-dot-dashed line). The jetlike outflow has a helical structure. 
The maximum jet velocity is mildly-relativistic (~ 0.3 c). However, it is clearly larger than 
the escape velocity. 
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Fig. 9. — Time variation of (a) the mass flux of the jet (solid line) and the accretion rate 
(dashed line), (b) the kinetic energy (solid line) and the Poynting flux (dashed line), and 
(c) the jet (solid line) and Alfven wave components (dashed line) of the Poynting flux at 
z/r s ^ 12 for the standard case A2. 
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Fig. 10. — Strong magnetic field case A3 (B = 0.07, v = 0.01). (a, b) Snapshots of density 
and plasma beta at t/r s = 120. The color scales show the values of the logarithm of density 
and plasma beta. The white curves represent the magnetic field lines. Arrows depict the 
poloidal velocities normalized by the light velocity. The central black region corresponds 
to the inner boundary. The distance and the time are given in units of and t$ = rs/c, 
respectively, (c, d, e, f, g, h) Various physical quantities on the surface of x/r s = 11 at 
tjr s = 120. The jetlike outflow is weaker and fainter (the maximum poloidal velocity of the 
jet is < 0.05 c ) than in the standard case A2. The magnetic twist is also weaker. 
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Fig. 11. — Weak magnetic field case Al (B = 0.03, v = 0.01). (a, b) Snapshots of density 
and plasma beta at t/r s = 170. (c, d, e, f, g, h) Various physical quantities on the surface 
of x/r s = 10 at t/rs = 170. The jetlike outflow is stronger (the maximum poloidal velocity 
of the jet is ~ 0.1 c ) than in the standard case A3. The magnetic twist is also larger. 
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Fig. 12. — Dependence of (a) the maximum vertical velocity of the jet (v z ), (b) the maximum 
toroidal velocity of the jet (v^), and (c) the maximum ratio of toroidal to poloidal magnetic 
field component (B^/Bp) of the jet (solid line) and of Alfven waves (dashed line) on the 
initial magnetic field strength (-Bo)- 
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Fig. 13. — Results of the fast rotation case B2 (B = 0.05, v = 0.015). (a, b) Density and 
plasma beta (= P gas / ' Pmag) at t/rs = 170. (c, d, e, f, g, h) Various physical quantities on the 
surface of x/r$ = 10 at t/rs = 170. The jetlike outflow is stronger and wider (the maximum 
poloidal velocity of the jet is ~ 0.1 c ) than in the standard case. 
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Fig. 14. — Results of the slow rotation case Bl (B = 0.05, v = 0.005). (a, b) Density and 
plasma beta at t/rs = 150. (c, d, e, f, g, h) Various physical quantities on the surface of 
x/r s = 10 at t/rs — 150. A jetlike outflow is not generated. 
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Fig. 15. — Dependence of (a) the maximum vertical velocity of the jet (v z ), (b) the maximum 
toroidal velocity of the jet (f^), and (c) the maximum ratio of toroidal to poloidal magnetic 
field component (B^/Bp) of the jet (solid line) and of Alfven wave (dashed line) on the initial 
rotational velocity (y rot ). 
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Table 1. Units and Normalization 



Phisical Quantity 


Description 


Normalization Unit 


rs 


Time 


r s /c 


r, z 


Length 


rs 


P 


Density 


Po 


P 


Pressure 


Poc 2 


V 


Velocity 


c 


B 


Magnetic field 


vW 2 



Note. — The unit of length is the Schwarzschild radius 
r s = 2GM Q /c 2 . The unit of density p is the initial density 
at (r,z) = (3r 5 ,0). 



Table 2. Models and Parameters 



Case 


B 


^0 


Eth 




E 

i-'mag 




E ro t 




Al 


0.03 


0.01 


2.36 x 10" 


-3 


6.04 x 10- 


-4 


5.36 x 10- 


-2 


A2 


0.05 


0.01 


2.36 x 10- 


-3 


1.68 x 10- 


-3 


5.36 x 10- 


-2 


A3 


0.07 


0.01 


2.36 x 10" 


-3 


3.28 x 10- 


-3 


5.36 x 10- 


-2 


Bl 


0.05 


0.005 


2.36 x 10- 


-3 


1.68 x 10- 


-3 


2.68 x 10- 


-2 


B2 


0.05 


0.015 


2.36 x 10" 


-3 


1.68 x 10- 


-3 


8.03 x 10- 


-2 


C 


0.05 


0.0 


2.36 x 10- 


-3 


1.68 x 10- 


-3 


0.00 




D 


0.0 


0.01 


2.36 x 10- 


-3 


0.00 




5.36 x 10- 


-2 



Note. - Case A2 (B = 0.05, v = 0.01) is considered to 
be the standard case in our simulations. Cases "A" differ from 
case A2 in initial magnetic field strength and cases "B" differ 
from A2 in the initial rotational velocity. In case "C" there is 
no rotation while in case "D" there is no magnetic field. 



